A possible explanation for Earth's climatic changes in the past 

few million years 

W. Wolfli * and W. Baltensperger t 
February 2, 2008 



Abstract 


4 


Planet Z: a 


transient member of the 


solar system 






The astronomical theory of Milankovitch relates the 


1.1 The orbit of planet Z and its close 




changes of Earth's past climate to variations in in- 




encounters with the Earth 




solation caused by oscillations of the orbital param- 




1.2 Solar radiation effects: the gas cloud] 


eters. However, this theory has problems to account 




1.3 Tidal effects on Z 





for some major observed phenomena of the past few 



million years. Here, we present an alternative expla- 5 The climate impact of the gas cloud oi 


nation for these phenomena. It is based on the idea planet Z 


that the solar system until quite recently contained 
an additional massive object of planetary size. This 


5.1 A gas cloud driven Greenhouse effect] 


5.2 The climate changes during the Ice 


object, called Z, is assumed to have moved on a Age epoch 





10 



11 

11 



It influ- 



enced Earth's climate through a gas cloud of evap- 6 Earth's past climate: a record of planet 


orated material. Calculations show that more than Z's history 






once during the last 3.2 Myr it even approached the 
Earth close enough to provoke a significant shift of 


3.1 Some more polar shifts 




3.2 On the origin and fate of planet Z . 



the geographic position of the poles. The last of 
these shifts terminated the Earth's Tee Age epoch. 



7 A new understanding of Earth's past 



14 



15 



15 



The origin and fate of Z is also discussed. 



climate 



Contents 



1 The problem of Earth's past climate 



changes 



2 Evidences for a polar shift event at the 



end of the last Ice Age period 



3 


The excitation of a polar shift during a 




close flyby of a massive object 




3.1 Earth's deformation in a strong grav- 




itational field gradient 




3.2 The fast polar 


wandering 




3.3 The impact of the polar shift event 




on the climate 





♦Institute for Particle Physics, ETHZ 
Honggcrbcrg, CH-8093 Zurich, Switzerland; e-mail: 
woelfli@particle.phys . ethz . ch 

^Centro Brasileiro de Pesquisas Fi'sicas, r.Dr.Xavier 
Sigaud,150, Urea, 22290 Rio de Janeiro, RJ, Brazil; e-mail: 
baltens@CBPF.br 



1 The problem of Earth's past 
climate changes 

Ancient climate is recorded in deep sea sediments, 
polar ice and some other natural archives mainly 
through changes in the oxygen isotope ratio, the 
S 18 0. In deep sea sediments the 5 18 0-signal reflects 
the combined effect of water temperature and of the 
amount of ice on Earth, whereas in polar ice it is a 
reliable proxy for the ambient local surface air tem- 
perature at the time of freezing. The study of pa- 
leontology and isotopic composition of hundreds of 
drilling cores extracted from these archives during 
the last few decades confirm that Earth's past cli- 
mate changed between colder and warmer phases, 
known as the glacial-interglacial cycles. Spectral 
analysis of this kind of proxy data covering the Late 
Pleistocene have yielded periodicities, which agree 
more or less with those of Earth's main orbital pa- 
rameters, i.e. precession (23 kyr), obliquity (41 kyr) 
and eccentricity (100 kyr), which appear at first 
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sight to support the old idea of Milankovich Q, 
that Earth's climate reflects the variation in inso- 
lation (solar heating) caused by the combined effect 
of the variation of these parameters ^ E|. How- 
ever, this astronomical theory of climate has sev- 
eral shortcomings. In particular, the main glacial- 
interglacial transitions occurred approximately ev- 
ery 100 kyr. This feature contradicts the astronom- 
ical calculations which show that the modulation of 
the insolation due to the eccentricity oscillation is 
negligibly small compared to that of the two other 
orbital parameters. Many hypotheses have been for- 
mulated to explain this "100 kyr problem", some- 
times involving additional astronomical parameters, 
such as the oscillation of the inclination of Earth's 
orbital plane, whose periodicity equals that of the 
eccentricity §, §. Most of the more recent at- 
tempts to understand this problem rely either on 
a non-linear response of the ice-sheet dynamics || , 
or of that of the coupled Atmosphere-Ocean sys- 
tem to the astronomical insolation forcing fit]. Al- 
though the results of some of these models compare 
well with the time spectra of the climate record, up 
to now none of them could reproduce the ampli- 
tude and phase of each glacial cycle. For instance, 
one of the most prominent intcrglacial periods of 
warmings observed around 400 kyr ago, known as 
the "isotope stage 11"- problem, happened at the 
time when the solar insolation was smaller than dur- 
ing the glacial times prior to and after this period. 
More recent measurements with improved time res- 
olutions, on deep sea sediment |8[ |t| and polar ice 
cores, GRIP @, GISP2 (HJ and Vostok jjJJ re- 
vealed some additional and most puzzling features 
of past climate. Fig. |l| shows the well known <5 18 0- 
record of the GRIP ice core covering the last 130 kyr 
of Earth's climate history, including the last inter- 
stadial, the Eem, the last glacial maximum around 
about 18-20 kyr BP, and its abrupt termination 
at about 11.5 kyr ago. Although some discrepan- 
cies between the GRIP and GISP2 data older than 
about 100 kyr BP have been noted pj| , it is evident, 
that Earth's climate changed rapidly during the last 
Ice Age period. These temperature peaks, called 
Dansgaard-Oeschger events, cannot be explained in 
terms of the Milankovitch effects, since they are too 
fast in comparison with the orbital parameters in- 
duced solar insolation variations. Even more puz- 
zling is the fast termination of the last Ice Age pe- 
riod, the Pleistocene-Holocene transition, shown in 
the lower part of Fig. [l| This transition ended the 
Pleistocene, and the following Holocene is a period 
of warmer temperature and generally more stable 
climate. It has been demonstrated that the climatic 



system in the Northern hemisphere is capable to 
switch rapidly from one given state to another |l3}| . 
But again, if for the same reason as mentioned above 
the Milankovitch effect cannot be the cause of such 
a switch, what else could it be? Is it possible that 
the luminosity of the sun varied in accordance with 
these observations? This is very unlikely, because 
prior to the Ice Age epoch, i.e. during the Late 
Pliocene until about 3.3 Myr ago, the climate was 
mild, stable for millions of years, and very similar to 
what it has been since the beginning of the Holocene 
l| . Why should the irradiance of the sun suddenly 
get noisier for a limited time? The Milankovitch 
effect and small changes of the solar irradiance are 
of relevance for Earth's climate. However, we agree 
with some of the critics pointing out that both pos- 
sibilities fail to explain many properties of the ob- 
served past changes of climate and that new ideas 
are required to solve the riddle of the Ice Age j| |(| . 

In what follows we present such a new idea. 
First, we show that the fast termination of the last 
Ice Age period, the Pleistocene-Holocene transition, 
can best be understood in terms of a rapid po- 
lar wandering, as has been suggested (and refused) 
more than 100 years ago. We then demonstrate that 
such a shift can occur during a close encounter be- 
tween an object of planetary size and the Earth. As- 
suming that this object, called Z, was on a highly 
eccentric orbit bound to the sun, we find that it 
could be responsible for Earth's climatic changes 
over the last few million years. Z must have dis- 
appeared from the solar system some time after the 
polar shift event, and we also consider the mecha- 
nisms responsible for its disappearance. 



2 Evidences for a polar shift 
event at the end of the last 
Ice Age period 

The Last Glacial Maximum (LGM) around 18-20 
kyr BP is characterized by the lowest known tem- 
perature in Earth's history, low precipitation rates, 
reduced sea level of at least 120 Meter |l4| , and huge 
polar ice shields extending over the northern part of 
America down to a latitude of almost 40° N during 
the last glacial maximum (Wisconsin stage), and 
in Europe to about 50° N, known as Weichsel and 
Wiirm stages in its northern part and in the Alps, 
respectively |]l5| . It is also well documented that 
at the same time large flocks of Mammoths, Woolly 
Rhinos and many other species of mammals roamed 
in north-eastern Siberia and north-western Alaska 
Il6l. Remains of these animals have been found on 
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Figure 1: Upper part: Ice core <5 18 record (GRIP) over 150 kyr 1 1 . This time span covers the last 
Interstadial, the Eemian, the last glacial period (ca. 90 - 11.5 kyrs) which is characterized by transient 
and fast temperature variations (21 Dansgaard-Oeschger events), and the Holocene reflecting the sudden 
improvement of Earth's climate at about 11.5 kyrs ago. The lower part is an expanded view of the same 
record for the last 20 kyrs. It shows the time structure of the Dansgaard-Oeschger event No.l and the 
fast Pleistocence-Holocene transition as indicated by the vertical dashed line. 
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the New Siberian Islands as far as 76° North in- 
dicating that food was available to support these 
largest mammals throughout the Late Pleistocene. 
It has also been noted, that at the same time this 
part of the world was covered by non-arctic forests 
up to about the same latitude, e.g. more than 1000 
km within present days' arctic cycle (66,66°). Even 
more puzzling, however, is the evidence that at the 
end of the last ice age, the so called Pleistocene- 
Holocene transition, these animals were not only 
instantaneously killed by some yet unknown cause, 
but also immediately deep frozen together with the 
surrounding mud p7| . The fact that the flesh of 
some of these animals was still eatable (for dogs) 
at the time of their discovery about 200 years ago 
suggests temperatures below freezing point through- 
out the Holocene in this part of the world. This is 
in contrast to the general warming up observed in 
many other places on Earth including both polar 
regions. The observation of simultaneous mass ex- 
tinctions of small and large animals in the western 
part of Alaska, in Western Europe and North Amer- 
and the distinct sea level rise leave no 
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doubt that the observed fast rise in temperature at 
the end of the last ice age period was accompanied 
by a global cataclysm and a rapid reorganization of 
the global climate. 

During the LGM the northern polar ice cap was 
not centered at the position of the (present day) 
north pole, but displaced by about 18-20 degree 
southwards on the American side of the globe Jl5| . 
This introduced the possibility of a rapid polar wan- 
dering, i.e. a change in orientation of the Earth with 
respect to its rotation axis. This idea was consid- 
ered in detail more than 100 years ago by the most 
eminent scientists of that time such as G.H. Dar- 
win, J.C. Maxwell, G.V. Schiaparelli, W. Thom- 
son and many others. They concluded that such 
a turn of the globe changing the geographical posi- 
tion of the North Pole to the amount required by 
the observations is in principle possible, but would 
require a force sufficient to overcome the stabiliz- 
ing effect of the equatorial bulge of our planet, 
which acts like the rim of a gyroscope. Since they 
were unable to find any terrestrial or extraterres- 
trial mechanism that could give rise to such a force, 
this idea was finally dismissed as impossible and, 
in fact, "not worth discussing" EjJ. Only more re- 
cently was this verdict seriously questioned when 
the evidence from rock magnetism studies unmis- 
takably suggested that the wandering of the ge- 
ographic poles was more rule than exception in 
Earth's history J2l], [22]] . However, this wandering 
is a slow process, probably connected to the conti- 



nental drifts first predicted by Wegener |23|], and 
therefore it cannot explain the rapid termination 
of the last ice age. In view of this dilemma, Hap- 
good Q suggested an alternative possibility based 
on the idea that the observed changes reflect the 
occasional sliding of the whole rigid crust over the 
viscous, plastic or possibly fluid inner layers of our 
Earth. He showed that the mass of the Antarctic ice 
cap is not symmetrically distributed about the rota- 
tion axis and believed that the resulting centrifugal 
effect is sufficiently large to initiate such a slippage, 
provided the ice cap can re-grow large enough after 
each slide. In his foreword to Hapgood's book, A. 
Einstein cautiously commented this idea as follows: 
"The only doubtful assumption is that the Earth's 
crust can be moved easily enough over the inner 
layers" . This is indeed the crucial point of this idea 
which eventually has not been accepted. 

Was this the end of the polar shift debate, al- 
though a fast polar shift is still the most straightfor- 
ward explanation for all observed effects connected 
with the fast termination of the last Ice Age period? 
In further investigating this question, we found that 
one particular excitation mechanism, namely the de- 
formation of the Earth in a transient external grav- 
itational field gradient, can produce the required ef- 
fect. 



3 The excitation of a polar 
shift during a close flyby of 
a massive object 

Any redistribution of Earth's mass is accompanied 
by a wobble and turn of the entire Earth relative 
to its rotation axis. This movement depends on 
both the initial mass redistribution and the rate 
at which the Earth's rotational bulge can readjust 
to the changing position of the rotation axis on 
the globe. Generally, the polar wandering resulting 
from Earth-bound causes are either too slow (drift 
of the continental plates), or too small (storm, ocean 
currents), to produce the rapid and large shift re- 
quired to explain the fast termination of the last ice 
age period. But what about cosmic events? There 
are two possibilities worthwhile to be considered: 
A collision of an object with the Earth, and a close 
flyby of an object. The first case has previously been 
discussed in the Literature It can be dismissed, 
because an impact, capable to significantly alter 
Earth's angular momentum or its moment of iner- 
tia, would have left one or several craters and other 
distinct geophysical and chemical signatures, which 
would have been discovered by geologists, bearing in 
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mind that the event we are discussing occurred only 
11.5 kyr ago. We therefore consider the possibility 
that the Earth had a close encounter with a massive 
object. We estimate first the minimum amount of 
deformation required to cause a polar wandering of 
about 18°, as is mentioned in chapter 2. We then 
show that such a shift can be induced by an object 
of planetary size during a flyby and that it is fast 
enough to explain the rapid termination of the last 
Ice Age epoch. 

3.1 Earth's deformation in a strong 
gravitational field gradient 

Prior to any disturbance the rotating Earth is to a 
large extent in a hydrostatic equilibrium. The mo- 
ment of inertia is a symmetric tensor S which con- 
nects the vectors of angular velocity uj and angular 
momentum D by the equation D = Hu;. For the 
undisturbed Earth the diagonalized tensor has the 
elements H[l, 1] = A = E[2, 2]=B = 8.010247- 10 37 
kg m 2 and 3[3, 3] = C = 8.036559 • 10 37 kg m 2 ||. 
The ratio C/A = 1.0032847 shows that the equa- 
torial bulge contributes about 3.3 %o to the main 
principal moment of inertia; the radius at the equa- 
tor exceeds that at the poles by about 21 km. This 
equatorial bulge stabilizes the position of the rota- 
tion axis along the figure axis C. During the flyby of 
a massive object, the globe is transiently deformed 
so that the new principal axis of the inertial tensor 
S — the figure axis — deviates from the rotation 
axis. These axes then precess around each other, 
whereby the angular momentum D remains fixed in 
space. For an observer on Earth the angular ve- 
locity vector uj moves around the figure axis. In 
co-ordinates fixed to the Earth this free motion is 
described by the Euler equation 



dE 



dt 



(1) 



where the square bracket represents the vector prod- 
uct. According to Eq. [l]the precession frequency of 
a rigid Earth would be wo(C — A) /A, with uio, the 
angular velocity of Earth's rotation. This yields an 
"Euler period" of 305 d. A minute precession of the 
rotation axis was observed by S.C. Chandler more 
than 100 years ago, but its period turned out to 
be 435 d Newcomb @ showed that this dis- 
crepancy can be understood by assuming that the 
Earth is not rigid but can be plastically deformed, 
which was quite a pioneering conclusion in those 
days. Since then, this so called Chandler's free mo- 
tion of the rotation axis has been studied in great 
detail, but always in view of small deviations only 



28 , 29J . The presently continuing Chandler motion 
is still open to discussion poj . 

A massive object passing the Earth at a close 
distance deforms it with a time dependent tidal 
force, which is the difference between the gravita- 
tional force at the given point and that acting at 
the Earth's center. When the distance R between 
the centers of the two masses is large compared to 
the radius Re of the Earth, the tidal force F(z) is 
parallel to the z-direction which points to the per- 
turbing mass Mz- It has the value 



F(z) = 2M z Gz/R 3 



(2) 



where G = 6.673 • 10 -11 m 3 kg _1 s -2 is the gravita- 
tional constant. On Earth's surface z = i?£COS7, 
where the angle 7 is the latitude. This tidal force in- 
creases as the third reciprocal power of the distance 
R between the center of the perturbing mass and 
that of the Earth. For example, if the moon would 
be 10 times closer to the Earth, then its tidal force 
would be 1000 times larger than it is at present. 

In first order approximation, this force tends to 
give the Earth an ellipsoidal shape. Such a deforma- 
tion is possible because Earth's interior is liquid or 
plastic except for the inner core and the peripheral 
crust [^5|. We estimate the deformation due to the 
force, Eq. |^, in a static approach by minimizing the 
energy as a function of the deformation in the exter- 
nal field. For this we assume that the redistribution 
of material at the surface of the Earth is given by 



H( 7 ) = # (cos(2 7 ) + 1/3) 



(3) 



where Hq is the amplitude of the tide. The energy of 
the material in the deformed state minus that prior 
to the deformation in the presence of the perturba- 
tion is then 



AE(H ) = ^pR% ( II ~, 



which is minimized with 



HqMzGR e 
R? 



H = R 2 E M z G/(2R 3 g) 



(4) 



(5) 



Here g = 9.8 m/s 2 is the gravitational accelera- 
tion at the surface of the Earth and p the specific 
weight of the displaced material(3500 kg/m 3 for the 
Earth's surface). Applying Eq. gto the Earth-Moon 
system we find that the peak-to-peak level differ- 
ence on Earth induced by our Moon is 2H$ = 0.35 
m which agrees quite well with the recently evalu- 
ated deformation of about 0.30 m [|lj. We use this 
formula to estimate the deformation produced by 
a flyby of a massive object, although we are aware 
that this simple approach probably overestimates 
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the deformation remaining immediately after such 
a flyby. In particular, if the response to the strong 
tidal force is too fast, then part of the deforma- 
tion will quickly disappear when the force ceases. 
We also expect that the actual deformation will not 
be a linear function of the gravitational force, and 
will depend on the relative motion of the perturbing 
mass to the surface of the rotating Earth. 

To estimate the minimum amplitude Hq of the 
tidal bulge required to cause a polar wandering of 
about 18°, as explained in chapter 2, we added to 
the equilibrium deformation of the Earth an addi- 
tional mass distribution of the form of Eq. || which 
was rotated by an angle a with respect to the orig- 
inal figure axis. The diagonalization of the new in- 
ertial tensors resulted in a new figure axis whose 
tilting angle relative to the former axis depends on 
Hq and a. Minimizing Hq as a function of a, we 
find that a tilting angle of 18° can be achieved with 
a bulge amplitude of Hq = 6.45 km and an angle a 
= 30°. This result allows us to estimate the mass of 
our ffyby-object or its distance of closest approach, 
by means of Eq. || which fixes Mz/R 3 - Assuming 
for instance that the distance of closest approach 
should be larger than the Roche limit, which for 
our Earth is about R t = 2.UR E , i.e. 15'600 km 
(assuming equal density of both objects), then the 
mass Mz must be larger than about 0.03 times the 
mass of the Earth (Me), a value which is somewhat 
smaller than that of Mercury. On the other hand, 
2-body dynamics tells us that the mass has to be 
much smaller than that of our Earth. Otherwise 
Earth's orbit, in particular its inclination and ec- 
centricity, would have been altered during such an 
event to an amount which would be incompatible 
with the present orbital parameters of the Earth. 
This consideration suggests an upper limit of about 
Mz < 0.2Me, and, according to Eq. |[ a maximum 
distance of closest approach of 29 300 km. Since our 
static approximation most probably overestimates 
the tidal effects, we suspect that the mass of Z was 
larger than 0.03 Me and that its closest approach 
during this event was even within the Roche limit. 

3.2 The fast polar wandering 

The sudden deformation of the Earth resulting from 
a short flyby event (the excitation time is less than 
10 h) defines the new figure axis around which the 
rotation axis immediately starts to precess, as de- 
scribed by Euler's Eq. |[ Since the figure axis it- 
self moves only by small amounts in co-ordinates 
fixed to the Earth, the rotation axis effectively spi- 
rals to the shifted position of the figure axis. To 
demonstrate this polar wandering, we solved Eq. [l] 



numerically. Since the Earth is not a rigid body 
we started the calculation with the inertial tensor, 
as defined by the initial deformation, and assumed 
that it relaxes towards the equilibrium tensor So 
oriented along the instantaneous rotation axis as 

cE = 5(f) -5p [3(f)] 

dt t [ ' 

Thus we approximated the relaxation behavior by 
a single relaxation time r, which includes both the 
capability of the equatorial bulge to follow the rota- 
tion axis and the disappearance of the deformation 
responsible for the polar wandering. Assuming a 
relaxation time of r = 500 d, we find that the ro- 
tation axis moves into the new position with a pre- 
cession period of about 400 d. The final destination 
is reached within about four precession periods or 
about 4.8 years (Fig. ^). 

This is certainly fast enough to explain the rapid 
termination of the last Ice Age period shown in the 
lower part of Fig. |l|. We have to add here, that the 
true relaxation behavior of the Earth in the presence 
of a large, but short-term deformation is still open 
to speculation. The theoretical estimates for this 
relaxation time vary between a few 10 d and 10 yr 
HH . We find that r in our example has to be larger 
than about 100 days to achieve a polar shift of the 
required magnitude. 

3.3 The impact of the polar shift 
event on the climate 

The lower part of Fig. [l] shows, that the termina- 
tion of the last ice age period is characterized by 
two events: A short term increase of temperature 
around about 14.5 kyr. ago, which was followed by 
a slow decrease back to the Ice Age level. The fol- 
lowing cold period, known as Younger Dryas (YD), 
lasted for about 800 years. At about 11.5 kyr ago 
the temperature suddenly rose again, but, and in 
contrast to the previous event, it continued to in- 
crease slowly until about 9000 years ago. From then 
on the surface temperature over Greenland essen- 
tially remained at that high value up to the present 
day. Can we understand why the first temperature 
increase was a transient episode only, whereas the 
second event had such a long lasting effect? To 
answer the first part of this question we need to 
consider an additional property of our flyby object, 
which will be discussed in the next chapter. It 
is sufficient to mention here that in this case the 
flyby distance of our object was simply too large 
to excite a significant polar wandering. The con- 
sequences of the second event can be explained in 
terms of the polar wandering discussed above. For 
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Figure 2: Left side: Geographic movement of the rotation axis induced by the fiyby of a massive object, 
assuming the deformation and relaxation parameters discussed in 3.1. and 3.2. The circles are latitudes 
with respect to the previous pole position, spaced 15°. Right side: Time dependence of the angle 8 
between the angular momentum (fixed in space) and the previous pole direction. 



this we have to remember (see chapter 2) that the 
shape and extension of the ice shield in the North- 
ern hemisphere prior to the shift suggest a pole po- 
sition at a longitude of about 90° W and a latitude 
of about 72° N, which is quite close to the present 
position of the North magnetic pole. Identifying 
this position with the starting point of our calcula- 
tion presented in Fig. |[ we see that the net effect 
of this polar wandering was a rotation of the globe 
by about 18°. As a consequence, the latitudinal 
pattern of the incident solar radiation immediately 
changed, so that, e.g. the low latitude part of the 
North American and European ice shields rapidly 
melted away, whereas in formerly ice free regions of 
the Arctic Ocean permanent ice was formed. Simi- 
larly, in the eastern part of the Antarctic the extent 
of the ice shelf was reduced, and, as we suspect, the 
formerly ice free western coast was covered to the 
extend observed today. After a transition time of 
about 2 000 years during which the polar ice sheets 
adjusted to the new pole position, the new state of 
the climate, called Holocene, was established. Com- 
pared to the previous climate state, the Holocene is 
significantly warmer, essentially, because the North 
pole moved away from the upper part of the North 
American continent into the present position in the 
Arctic Ocean. This resulted in an increased heat 
exchange between ice and water, which significantly 
reduced the extent of the polar ice cap. The albedo 
decreased accordingly. This reduction was only par- 
tially compensated by the corresponding changes in 
the southern hemisphere. According to this sce- 
nario, the Holocene is a new stable epoch and not 
just another Interstadial, as is generally believed, 



because the new pole position is stable and Z no 
longer influences the Earth. 

4 Planet Z: a transient mem- 
ber of the solar system 

We have seen that an object of mass Mz between 
0.03 Me and 0.2 Me can excite a polar shift which 
is both fast and large enough to explain the termi- 
nation of the last Ice Age, provided it passed the 
Earth at a distance of less than about 30 000 km. 
Since it is extremely unlikely that an object on a 
hyperbolic trajectory would accidentally approach 
the Earth so close during one singular fiyby, we con- 
sider the possibility that it was moving on a bound 
orbit extending beyond that of the Earth. The prob- 
lem with this assumption is that this object, Z, no 
longer exists — otherwise we would observe it — so 
that the occurrence of the proposed close encounter 
at a specific time cannot be verified by backward 
calculations. Therefore, we shall consider the long 
term behaviour of Z on the basis of assumed orbital 
parameters in order to find out whether close ap- 
proaches to the Earth occur as often as expected. 
We shall see that this "trial and error" approach 
not only gives a positive answer, but even provides 
an explanation for the frequent fast climatic changes 
observed prior to the pole shift event (see chapter 
5). 
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4.1 The orbit of planet Z and its close 
encounters with the Earth 

There are constraints on the initial orbital param- 
eters of Z. Its aphelion has to be larger than that 
of the Earth, otherwise it could not approach the 
Earth close enough. We shall see further below that 
at each perihelion passage a surface layer of Z has to 
evaporate at sufficiently high temperatures so that 
atoms reach the escape velocity. This requires an 
extremely eccentric orbit with a -C b, where a is the 
perihelion and b the aphelion distance. Thus the 
major semi-axis A = |(a + b) f» ib must be larger 
than 0.5 AU, and the eccentricity e = (1 — a/A) close 
to one. Expecting that the encounter frequency 
rapidly decreases with increasing A, we restrict the 
investigations, to the range 0.5 < A < 2 AU. There 
are no restrictions on the inclination i of Z's orbit, so 
that the full range 0° < i < 180° has to be explored. 

With this in mind, Nufer et al. (3j| calculated 
the movement of Z in the gravitational held of the 
sun, the Earth, the Moon and all other planets, ex- 
cept Uranus, Neptune and Pluto, and evaluated by 
means of a Pascal program its closest approaches to 
the Earth for each orbital period. Details of this 
program, its performance, and the relevant results 
shall be published in f$3j. In the simulations Z 
was always added to the solar system at the year 
J2000.0, ignoring the fact that at this time it no 
longer existed. So far calculations have been per- 
formed, forward as well as backwards in time, for 
two different masses (Mz — 0.11Mb equal to that 
of Mars, and Mz = 0.06Mb) and seven different 
sets of initial orbital parameters (A — 0.978 and 2 
AU, six different inclinations in the range i = 0°- 
180° and two eccentricities e = 0,973 and 0.987) 
The time ranges varied between 150 and 750 kyr. 
As expected, we found that encounters are more fre- 
quent for A = 0.978 than for A — 2. Their number 
also increases with increasing eccentricity e, and is 
largest for small inclinations i < 6°. The inclina- 
tion is defined here as the angle between Z's orbital 
plane and the invariable plane which is perpendic- 
ular to the orbital angular momentum of the solar 
system. The results obtained with the initial condi- 
tions A = 0.978, e = 0.973, i = 0° and 180°, respec- 
tively, are too similar to decide whether Z was on 
a pro- or a retrograde orbit. Since the computing 
speed strongly depends on the precision assumed for 
the numerical integration procedure, we restricted it 
to a moderate value of 10 -13 per integration step, 
which was sufficient for most considerations. With 
this restriction, about 15 h were needed to cover 10 
kyr using a 133 MHz PC. 

The result of one of these calculations is dis- 



played in Fig. |3|. It shows the closest approaches 
of Z (Mz = 0.11Mb) to the Earth over a time pe- 
riod of 750 kyr. The initial parameters at J2000.0 
were A = 0.978 AU, e = 0.973 and i = 0°. Plot- 
ted are all encounters with distances of less than 
0.02 AU = 3 • 10 6 km. The irregular clustering of 
these approaches is the result of the osculating or- 
bital parameters of Z and the Earth. An expanded 
view of all closest approaches with distances less 
than twice that of the Moon over 750 kyr is given 
in Fig. |]. To illustrate the mass dependence we 
also show the pattern obtained for an object with 
Mz = 0.06Mb;, but otherwise identical orbital pa- 
rameters. The two horizontal dashed lines in both 
figures indicate Moon's distance (384000 km), and 
the critical distance below which significant polar 
shift events have to be expected (30 000 km). In 
the first case we note a closest approach of about 
59 000 km, i.e. no pole shift event, but a distinct 
one in the second case (23 500 km). From these and 
the other trials we conclude that an object moving 
on a highly eccentric orbit approaches the Earth to 
polar shift distance about once in 1 Million years, 
which is compatible with the duration of the ice age 
epoch (see chapter 6.1). It is meaningless that the 
calculated polar shift event does not occur at the ex- 
pected time in view of the fact that the exact mass 
and orbit of Z is unknown. In both cases we also 
observe about 6 fly by 's per 100 kyr below 300000 
km. In all these events the gravitational interaction 
is still strong enough to excite strong earthquakes 
and volcanic activities (see chapter 7). All these 
close flyby's also influence Moon's orbit. We found, 
however, that the resulting disturbances over the in- 
vestigated time ranges are averaged to the amount 
compatible with the present orbit p3. The fre- 
quency of the flyby's increases with their distance. 
For distances larger than twice that of the Moon 
they are harmless regarding tidal effects induced on 
Earth, but they can still influence Earth's climate. 
To understand this additional effect, which is fur- 
ther discussed in chapter 5, we next examine what 
happened to Z on the proposed highly eccentric or- 
bit. For convenience, we assume that it was a Mars 
sized object with M z = 0.11Mb = 6.58 • 10 23 kg, 
and radius Rz = 3.39 • 10 6 m. 



4.2 Solar radiation effects: the gas 
cloud 

An object moving on a highly eccentric orbit is ex- 
posed to intense solar radiation. The total amount 
of energy Ez s received by Z per unit of surface and 
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Figure 3: Upper panel: Closest approaches of planet Z (Mz = 0.11Mb) to the Earth over 750 kyr and 
to distances of less than 0.02 AU = 3 • 10 6 km as indicated by the endpoint of each vertical line. Lower 
panel: Expanded view of the irregular clustering of these approaches over 150 kyr. 
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Figure 4: Closest approaches over 750 kyr below twice the Moon-Earth distance. The pattern depends 
on the selected orbital parameters and the mass of the object. The left side shows the result obtained 
for a mars- like object (Mz = 0.11Mb), and the right side that for an object with Mz = 0.06Me but 
otherwise identical start parameters. The two horizontal dashed lines in both figures indicate Moon's 
distance (384 000 km), and the critical distance below which significant polar shifts have to be expected 
(30 000 km), respectively. The open (filled) circles mark the encounters of the objects on their way towards 
(from) the sun. 
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orbital period is 



E Zs = Ps 



(It 



AirR 2 



Mz 
2L 



-P, 



(7) 



where P s = 3.73-10 26 J/s is the total power radiated 
by the sun. The orbital angular momentum L = 
M Z R 2 ^§ {M z - mass of Z, R its distance to the 
sun, and ip the angle co-ordinate) is determined by 
the major semi- axis A and the eccentricity e 



L = Mz^jM s GA(l-e 2 ) 



(8) 



G is the gravitation constant, and Ms the mass of 
the sun. Therefore 



Ezs — 



1 



2 y/M s GA{l - e 2 ) 



(9) 



Using the same orbital parameters as in Fig. || 
{A = 1 AU, e = 0.973), this formula yields an en- 
ergy per unit surface of 1.67 • 10 11 J/m 2 per orbit. 
During each perihelion passage the power density 
reaches a maximum of about 2.83 MW/m 2 . This 
is sufficient to evaporate any type of rock on the 
exposed side and also to heat up the resulting gas 
so that its light components, mostly atoms, even 
reach the escape velocity v e , which is 5020 m/s for 
a mars-like object. If all the radiation energy were 
converted into kinetic energy, then the total mass 
per m 2 being accelerated to escape velocity would 
amount to 

m g = 2^% (10) 
v 2 

or m g = 13 300 kg/m 2 . This corresponds to an abla- 
tion rate of almost 5 m per orbital period, assuming 
that Z was composed of silicates or basalts with typ- 
ical densities of about 2700 kg/m 3 . Since the con- 
version efficiency was hardly 100% , this represents 
an upper limit for the total mass of evaporated ma- 
terial. It amounts to M g = TrR 2 z m g — 4.8 • 10 5 Gt 
per orbital period, indicating that Z was surrounded 
by a rapidly expanding gas cloud, with Oxigen and 
Silicon (atomic weights 16 and 28, respectively) as 
its main components. Its radius at the distance of 
Earth's orbit was R g = T g v g , where T g is the travel- 
ling time from the sun to Earth's orbit, and v g the 
expansion velocity of the gas cloud. The tempera- 
ture T of the atmosphere of Z from which the fastest 
atoms escape equals the temperature of the escaped 
gas cloud. A measure of the radial expansion veloc- 
ity is ^J<v 2 > = 

constant, and m a the atomic mass of the gas com- 
pound. Assuming a temperature of T = 2000 K 
during the perihelion passage, this yields a veloc- 
ity of 1000 m/s for Oxygen atoms. In r g = 34 



hsIL^ w ith ks the Boltzmann 



days, the travelling time of Z from the perihelion 
to Earth's orbit, the gas cloud expanded to a ra- 
dius of 3 Mio. km. The cloud remained spherical 
during this time, provided the constituents of the 
cloud have no resonant transitions in the solar spec- 
trum. In this case, the Rayleigh scattering cross 
section is of the order of the square of the classical 
electron radius, so that less than one photon per 
atom was scattered during the travelling time. In 
this case the light pressure effect is negligibly small, 
and the scattered light is dominantly blue in colour. 
This behaviour differs from that of an opaque grain 
which is accelerated inversely proportional to its ra- 
dius by the light pressure. A particle of 1 micron is 
displaced by many Mio. km during 34 d. 

It is interesting to compare these considerations 
with observations: For instance, the coma of the 
comet observed in 1811 AD had a diameter of 1.8 • 
10 6 km J34|, comparable to the estimate obtained 
for Z. Its tail had a length of 179 • 10 6 km and was 
"white" in colour suggesting that it contained dust 
grains of various sizes and molecules with strong 
transition in the solar spectrum. The appearance of 
Z was different. Its nucleus was much larger than 
those of any presently known ordinary comets and it 
was surrounded by a huge but essentially spherical 
gas cloud with bluish colour. 

4.3 Tidal effects on Z 

Z was also exposed to strong time dependent gravi- 
tational forces during each perihelion passage. They 
induced a tidal wave, whose kinetic energy was con- 
verted into heat so that its temperature increased 
stepwise. To get an idea of this heating effect, we 
estimate first the amplitude of this tidal wave. For 
this we use Eq. [| but this time with M z for the so- 
lar mass, and Rz and gz for the radius and surface 
acceleration of Z {Mz, Rz, and gz equal to that of 
Mars). For a perihelion distance R — a = 4 • 10 9 m, 
as determined by the orbital parameters A = 1 AU 
and e = 0.973, we find a value of Hq = 3230 m. The 
lowering of the energy through the deformation in 
the tidal field, given by Eqs. 5 and 4, expressed in 
variables for Z, is 



2 R%GM% 

AE = z R s 

15 i? 6 



(11) 



where Ms is the solar mass, and R the distance 
to the sun. During an orbit of Z twice this energy 
is transformed into heat, i.e. 2\AE Z t\ = 1 • 10 25 
Joule for R = a. Tidal forces act on the sun 
as well. The corresponding values are 0.6 m and 
— 1 • 10 23 Joule, respectively, smaller than the tidal 
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effect on planet Z. Comparing AEzt with the in- 
flux of solar radiation which according to Eq. ^| is 
Ezs^Iiz = 6-05 • 10 24 Joule; hence the values are of 
the same order. Radiative energy primarily heated 
up the surface, while the tidal deformation extended 
into the interior of planet Z. Assuming that all of 
this energy was converted into heat, then the tem- 
perature increase was AT = 2 ^ E A z / z = 0.02 K per 
orbit, assuming a specific heat of C p = 650 J/kg 
K. After about one Million orbits, approximately 
equal to Earth years, the interior was melted and, 
according to the Stefan-Boltzman T 4 -law, tidal en- 
ergy alone could maintain a temperature of about 
450 K on a black surface. This implies that the 
surface was only molten during each perihelion pas- 
sage, whereas on its journey towards the aphelion 
a thin solid crust was formed on the surface. Once 
liquefied, Z easily reached its optimum shape during 
each perihelion passage. This redistribution of mat- 
ter was a turbulent process and occurred twice dur- 
ing each passage. For smaller perihelion distances 
tidal heating would dominate because of its a de- 
pendence per orbit as compared to the a -1 / 2 depen- 
dence of the radiation energy given by Eq. ^. 

Tidal losses imply that the deformation lags be- 
hind the exciting tidal field, so that the forces on 
Z are not symmetric on incoming and outgoing po- 
sitions. As a consequence, both the orbital energy 
Eq and the orbital angular momentum L decrease. 
With decreasing Eq the major semi-axis A of the or- 
bit is reduced. The effect of the tidal energy dissipa- 
tion can be estimated by comparing the maximum 
tidal loss 2AEzt with the orbital energy determined 
by E a = - MzMs£ = -2.65 • 10 32 Joule. The ratio 
7nr§ — = 4.4 • 10 7 is a measure for the number of 
orbits required for a substantial change of A. The 
influence on the eccentricity e is less evident. The 
relation e 2 = 1 + m^m^g 2 snows tnat a decrease in 
Eq lowers e, while a decrease in L increases e, since 
Eq is negative. In view of our restricted knowledge 
of Z, we are not able to state which of the two effects 
prevailed. 

5 The climate impact of the 
gas cloud of planet Z 

The upper panel of Fig. || indicates that Z frequently 
approached the Earth to distances within its sur- 
rounding gas cloud, which had a radius of about 
3 Mio. km, according to our estimate. We shall 
show that the incidence of this gas on Earth's at- 
mosphere resulted in a Greenhouse effect which was 
large enough to enhance the global temperature. 



Using a simple model for this Greenhouse heating 
we further demonstrate that the frequency struc- 
ture of this effect compares surprisingly well with 
that of Dansgaard-Oeschger events observed during 
the Late Pleistocene (Fig. Q). 

5.1 A gas cloud driven Greenhouse 
effect 

Each time the Earth passed through the gas cloud 
of Z, a fraction of this material was captured by 
Earth's upper atmosphere. For a spherical gas 
cloud, and constant gas density, the amount of gas 
m c captured during one crossing is proportional to 
the ratio of the volume covered by the Earth dur- 
ing its passage, V c < irR 2 E 2R g: and the volume of 
the gas cloud V g — 47ri? 3 /3. Neglecting the gravi- 
tational focusing, which enlarges the capture radius 
R c over Earth radius Re by at most 10 % , the total 
mass of captured particles is determined by: 

m c < 1.5m g R 2 E /R 2 g (12) 

With the estimated mass for the gas cloud m g = 4.8- 
10 5 Gt, its radius R g = 3T0 6 km, and Earth's radius 
Re = 6.4 • 10 3 km, this yields an amount of m c < 
3 Gt per encounter. The atoms of this cloud hit 
the upper atmosphere of the Earth with a relative 
velocity of about v r = 42 km/s, assuming a crossing 
angle of 90° and equal velocities of about 30 km/s 
for the Earth and Z. The particles are stopped at an 
altitude of about 150-200 km above Earth's surface. 
Including the acceleration in the gravitational field 
of the Earth (escape velocity v es = 11 km/s), we 
find that the total energy deposited at this altitude 
amounts to E g — m c (v 2 + v 2 s )/2 = 3 • 10 21 J. Since 
the passage through the cloud lasted up to 37 hours, 
this energy was deposited into a thin layer of the 
uppermost atmosphere of the whole globe. 

The impact of this cloud on Earth's atmosphere 
was manifold considering the fact that each Oxy- 
gen atom had a kinetic energy of about 160 eV. 
This energy was converted into radiation and used 
to ionise and dissociate the main constituents of the 
atmosphere, the Nitrogen and Oxygen molecules. 
The resulting ions and atoms initiated a complex 
chain of chemical reactions in the upper atmosphere 
during which the well known infrared-active com- 
pounds such as Ozone (0 3 ) and Nitrogen oxydes 
(NO x ) were formed. By means of the following sim- 
ple arguments we can get some idea of the total 
amount of infrared active molecules produced dur- 
ing such an event. Assuming that 30 % of the par- 
ticle energy, or 50 eV, is consumed by dissociation, 
then about 20 N2 or O2 molecules are decomposed 
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per incoming particle. The resulting 40 excited and 
ionised atoms interact with the molecules of the at- 
mosphere. Assuming that 50 % of these atoms end 
up in an infrared active molecule, then about 20 
of these would be produced per incoming particle. 
Since their mean molecular weight is about twice 
the atomic weight of Oxygen, we expect that up to 
20 ■ 2 • 3 Gt = 120 Gt of Greenhouse gases, mostly 
Nitrogen compounds and Ozone, could be produced 
during a flyby. True, this is less than the atmo- 
spheric inventory of CO2 which was of the order of 
430 Gt during the last glacial maximum. However, 
this deficit is more than compensated by the fact, 
that the radiative forcing of most of the infrared ac- 
tive molecules of interest here is 30-70 times larger 
than that of CO2 for equal concentrations. We 
conclude that the estimated amount of Greenhouse 
gases produced is sufficient to transiently enhance 
the mean global temperature to a level higher than 
what is predicted for a doubling of the CO2 con- 
centration [|5). The role of CO2 and CH 4 during 
such an event is open to discussion. According to 
the measurements made on air trapped in the po- 
lar ice cores, their concentrations essentially seems 
to follow the observed global temperature changes 
f3(| , suggesting that they play a passive role only. 

5.2 The climate changes during the 
Ice Age epoch 

In chapter 3 we have explained why the climate was 
generally colder during the Late Pleistocene than in 
the following Holocene. Here, we show that the sud- 
den increases in temperature of variable length and 
frequency, the Dansgaard-Oeschger events shown in 
Fig. |l|, can be explained in terms of the Green- 
house heating effect discussed above. For this we 
assume that the global temperature increase AT is 
proportional to the Greenhouse gas concentration, 
and that this value depends linearly on the num- 
ber of particles captured during each encounter. We 
further assume that the suddenly enhanced gas con- 
centrations decrease exponentially in time after each 
encounter. In this case, the time dependent temper- 
ature variations are given by 

AT(t)=T ^,/l-^|e-^ i (13) 
j V a 

were Dj is the distance of closest approach of the 
flyby object at time tj, R g the radius of the gas 
cloud, and r r an averaged residence time of the rel- 
evant Greenhouse gases. To is the initial tempera- 
ture increase for Dj <C R gi which depends on the 
amount of material captured from the gas cloud, the 



efficiency of conversion into Greenhouse gases, and 
the radiative forcing of these. The value of Tq is 
open to discussion. For the reasons explained above 
we are optimistic that Tq is large enough to account 
for the AT observed in the polar ice cores. We also 
note that the residence times of the various Green- 
house gases of interest here are not known. This 
crude model ignores long term reactions of Earth's 
climate system to the sudden changes in the radia- 
tion budget. 

Applying Eq. [[3] to the closest approach spec- 
trum given in Fig. || yields the Greenhouse heat- 
ing signal shown in Fig. ^[ Outgoing and incoming 
passages were equally weighted. Actually, incoming 
gas clouds are more expanded so that a single pas- 
sage contributes less. On the other hand they also 
contribute for larger distances than those consid- 
ered in Fig. [5| This calculation was performed with 
a time constant of r r = 200 yr which we consider 
an uppermost limit for the residence time of the 
Greenhouse gases of interest here. Again, a one-to- 
one correspondence between predicted and observed 
temperature spikes and clusters of spikes cannot be 
expected for the reason given in chapter 4. Nev- 
ertheless, a comparison of the calculated spectrum 
with the 5 18 polar ice data displayed in Fig. [I] re- 
veals striking similarities in shape, frequency and 
irregular clustering. E.g. the spectrum in Fig. ^ 
(lower panel) exhibits 17-26 distinct temperature 
excursions per 100 kyr, values, which compare well 
with the 21 Dansgaard-Oeschger events observed in 
the polar ice data (Fig. |l|) and some other palaeo- 
climate archives during the last glacial period, i.e. 
90-11.5 kyr ago §. 

In order to determine the frequency spectrum of 
the Greenhouse signal presented in Fig. ||, we ap- 
plied the Fast Fourier Transformation (FFT) algo- 
rithm to 2 19 data points of this spectrum, i.e. to 
the time range from +150 to -374 kyr. The result- 
ing power spectral density as a function of the re- 
ciprocal value of the frequency (the periodicity) is 
shown in Fig. ^|. Distinct maxima can be seen at 
6, 19, 23, and ca. 50, 90 and 175 kyr. These pe- 
riodicities are the result of the interplay between 
the complex movement of Z's and Earth's orbital 
planes and should not be confused with similar val- 
ues delivered by the Milankovitch theory. Here, the 
decisive parameters are the oscillation of both in- 
clinations (about 7 kyr, and 100 kyr, respectively) 
and the movements of both perihelia (33). In the 
high frequency domain these values compare quite 
well with those obtained from comparable analysis 
of the GRIP data displayed in Fig. [l|. The peak 
observed at 90 kyr is close to the famous 100 kyr 
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kyr since J20OQ.0 

Figure 5: Upper panel: Greenhouse spectrum for the distance of closest approach pattern displayed in 
Fig. |^. The vertical scale is in arbitrary units. Lower panel: Expanded view of this spectrum over 150 
kyr shows the similarity between the calculated structure and that of the Dansgaard-Oeschger events 
displayed in Fig. ||. 




Figure 6: Power spectral density (PSD) of the Greenhouse spectrum presented in Fig. ||. Visible are three 
distinct peaks at 6, 19, 23 kyr, and three broad maxima at ca. 50, 90 and 175 kyr. 
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observed in the deep see sediment data. This sug- 
gest that the Interstadials might be the result of 
an enhanced clustering of the Dansgaard-Oeschger 
events. It should be noted, however, that such a 
comparison is problematic for two reasons: First of 
all, the results of our spectral analysis depend on 
the selected orbital parameters and in particular on 
the mass of Z. Therefore, we can only say that the 
evaluated periodicities of the postulated Greenhouse 
effect are not in contradiction with observation. Sec- 
ond, any comparison with deep sea sediment data is 
questionable, because the time scale of these mea- 
surements is usually tuned to the 100 kyr cycle of 
Earth's orbital eccentricity. 

6 Earth's past climate: a 
record of planet Z's history 

In chapter 1 we pointed to the strange fact that 
Earth's climate was remarkably stable prior to the 
Ice Age epoch. Actually, the <5 18 record of benthic 
foraminifcras from the Ocean Drilling program, site 
659, shown in Fig. 0, suggests three different cli- 
mate regimes for the last 5.0 Myr B7j. The first 
phase is characterised by a very stable climate with 
small polar ice shields, and, in comparison with the 
following regimes, high temperature if it is assumed 
that ice extension and mean global temperature are 
correlated. This climate state is similar to that of 
the Holocene. A sudden drop in the global mean 
temperature terminated this state at about 3.2 Myr 
ago. In the following period the variability of the 
climate began to increase until about 1.0 Myr ago. 
At this time the proxy dates indicate a further drop 
of the mean temperature to the lowest known level 
in Earth history. This marks the beginning of the 
Ice Age epoch discussed in the previous chapter. It 
is obvious, that the properties of these three, or ac- 
tually four different climate epochs, including the 
Holocene, cannot be explained in terms of the orig- 
inal Milankovitch theory, unless changes of solar ir- 
radiance or Earth's orbital parameters according to 
observations are postulated. However, the introduc- 
tion of such changes is hard to justify. 

6.1 Some more polar shifts 

In view of what has been derived so far the spec- 
trum in Fig. ^ can be interpreted in a new way. 
Prior to about 3.3 Myr ago planet Z was inactive. 
If it existed, it was not yet on a highly eccentric or- 
bit, cold, and without gas cloud. The climate was 
stable showing only the moderate changes which ac- 
tually reflect the undisturbed Milankovitch effect. 



In agreement with this theory no 100 kyr-cycle ap- 
pears up to this time || . The fact that the size of the 
continental ice caps during this period was compa- 
rable to the present day value indicate that Earth's 
North pole was then located either in the North- 
ern part of the Atlantic, in the Pacific Ocean, or at 
about the same position as it is today. It should be 
possible to evaluate its actual position using infor- 
mation on the climate in various regions during the 
Pliocene, and/or from geomagnetic field measure- 
ments. Around 3.2 Myr ago we assume that planet Z 
shifted Earth's North pole towards the North Amer- 
ican continent in the way described in chapter 3. As 
a consequence, the continental ice shields started to 
grow and with it Earth's albedo, which reduced the 
global mean temperature. Of course, this implies 
that Z was on an eccentric Earth crossing orbit not 
later than at this time. Tidal heating consumed or- 
bital energy so that its major semi-axis was slowly 
reduced. As mentioned in chapter 4, we are unable 
to decide in which way its eccentricity changed. The 
growing amplitudes of the climatic changes between 
about 3.2 and 1.0 Myr suggest a slowly growing gas 
cloud, which indicates an increasing orbital eccen- 
tricity. The tidal effects increased as well, which 
accelerated the heating up of Z. About 1.0 Myr 
ago, a second polar shift in direction toward the 
North American continent to the position assumed 
in chapter 3 further lowered the global temperature. 
The Ice Age epoch began. At this time the gas 
cloud was developed to its full size and the resulting 
Greenhouse heating created the temperature varia- 
tions described in the previous chapter. A third and 
final polar shift terminated Earth's Ice Age epoch 
11.5 kyr ago, as discussed in chapter 3. The polar 
ice caps readjusted to the new position and, because 
of the diminishing size of the continental ice sheets, 
the global temperature increased during the early 
Holocene to its present value. 

The suggestion that Earth's pole was shifted at 
least three times during the last 3.2 Myr is compat- 
ible with the orbital calculations presented in chap- 
ter 4. They show that Z approached the Earth to 
polar shift distances about once in a million years, 
so that more than one shift event may be expected 
during the time range considered here. The direc- 
tion of the polar shift depends on the geographic 
position of the closest approach. Antipode events 
excite the same directional shift, while those at posi- 
tions differing by 180° in longitude only, or positions 
with the same longitude but equal northern, respec- 
tively southern latitudes lead to opposite motions. 
The earlier shifts occured with encounters outside 
the Roche limit, which are more likely and, in view 
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Age (Myr BP) 

Figure 7: Benthic foraminifera 5 1S record from site 659 (18° 05'N, 21° 02'), according to Ref.§7|- This 
spectrum is a proxy for the extension and variation of the polar ice shields over the last 5.0 My. Increasing 
S 18 corresponds to decreasing ice caps and to warmer climate. The late Pliocene (before 3.2 Myr) is 
characterized by small ice caps and remarkably stable climatic conditions comparable to those of the 
Holocene. The corresponding value is shown by the horizontal dashed line. 



of the larger mass of Z at early times, sufficient to 
produce the required deformation. 

6.2 On the origin and fate of planet 
Z 

At this point it is tempting to speculate on the ori- 
gin of Z. We see three possibilities: First, Z was 
originally on a hyperbolic trajectory but, during 
a close encounter with Jupiter, transferred to the 
required highly eccentric bound orbit. The major 
semi-axis A of such an orbit would be larger than 
half of Jupiter's orbital radius, i.e. A > 2.6 AU, 
because of the presumably large initial velocity at 
the inflection point. Such a large A is not compati- 
ble with our orbital calculations in chapter 4, which 
suggest A « 1 AU. Second, Z was originally a moon 
of Jupiter which was ejected due to the coupling 
with the other moons of Jupiter. Such an emission 
could result in a rather small initial velocity and 
accordingly to a highly eccentric orbit. Again, this 
is unlikely for the same reason as before. The last 
possibility relates to the Asteroid belt problem. It 
has been argued that the asteroids circling the sun 
at a mean distance of about 2.7 AU could be the re- 
mains of a collision between two larger bodies, one 
originally bound to the sun, the other presumably 
entering the solar system on a hyperbolic trajec- 
tory. Assuming that Planet Z was a reaction prod- 
uct of such a collision emerging with low velocity, 



then the major semi-axis of its orbit would be not 
much larger than about A = 1.35 AU at the begin- 
ning of its journey around the sun. Considering the 
time required to reduce the orbital energy and to in- 
crease its eccentricity to the final values estimated 
in chapter 4, it is the most likely scenario. 

It is open to discussion at what time Z disap- 
peared. One possibility is that in the last close en- 
counter it lost its angular momentum and plunged 
into the Sun. However, this requires an excentric- 
ity of not less than 0.98, and it supposes that the 
angular momentum received is opposed to the ex- 
isting one. Therefore it seems more likely that the 
last encounter with the Earth was within the Roche 
limit, so that Z split into several parts. As a result 
of the reduced escape velocities of the debris, these 
parts lost material at an accelerated rate and were 
completely dissolved during the Holocene. The gas 
and dust clouds were expelled from the planetary 
system by radiation pressure. 

7 A new understanding of 
Earth's past climate 

The previous chapters showed that the observed 
complex structure of Earth's climate changes dur- 
ing the last few Million years can be explained in 
terms of a single object of planetary size, which was 
bound to the sun on a highly eccentric orbit for a 
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limited time. This planet, called Z, was surrounded 
by a gas cloud and interacted gravitationally and 
through its gas cloud with the Earth and its atmo- 
sphere with a frequency determined by the pattern 
of its closest approaches to the Earth. These effects 
influenced Earth's climate in two different ways. 

Each polar shift event started with a huge flood, 
exceptionally strong winds, earthquakes and vol- 
canic activities and ended with a displaced pole po- 
sition after a couple of years. Following this shift 
the polar ice caps adjusted to the new pole posi- 
tion which changed Earth's reflectivity, the albedo, 
and accordingly the global temperature. The new 
pole position determined whether this temperature 
was higher or lower after the shift. In agreement 
with observation we found that this happened about 
three times during the last 3.2 Myr. Two times the 
temperature dropped, suggesting that each time the 
pole was shifted closer towards the North American 
continent, so that the size of the continental part of 
the polar ice cap increased. During the last event, 
11.5 kyr ago, it was moved into the arctic sea. At 
this position a net reduction of the polar ice caps 
prevailed. As a result of the reduced albedo, the 
global temperature increased to the presently known 
value. We also found that Z passed the Earth about 
6 times in 100 kyr at distances less than 300 000 km. 
All of these fly by 's were close enough to excite large 
tides followed by earth quakes and volcanic activ- 
ities. Recently, it has been suggested that the so 
called Heinrich events observed in the North At- 
lantic sediments during the last 100 kyrs — a rare, 
but distinct ice rafting phenomenon — could be the 
result of unusual, strong earthquakes ||| |39|| . It 
is interesting to note that the observed 6 Hcinrich- 
events per 100 kyr compare well with the predicted 
frequency for sub-lunar flybys. 

The interaction between the gas cloud of Z and 
Earth's atmosphere was quite spectacular as well, 
since most of the kinetic energy of the particles 
stopped in the upper atmosphere was converted 
into radiation. Probably, this radiation was in- 
tense enough to even ignite extended fires on Earth's 
surface. After a couple of days the fires were 
extinct either by the flood, by heavy rain or by 
snow enriched with nitrous acid. When the atmo- 
sphere was cleared of the light absorbing compo- 
nents, the radiative forcing of the produced Green- 
house gases began to prevail, so that globally the 
temperature rapidly increased within a few years 
to the maximum value determined by the amount 
of Greenhouse gases produced during the encounter. 
The Greenhouse effect was a transient phenomenon, 
whose impact was limited by the amount of the Ni- 



trogen compounds produced per encounter and the 
residence time of these gases. It may be responsible 
for the Dansgaard-Oeschger events observed during 
the Late Pleistocene. If so, then the Interstadials 
may be explained in terms of the quasi-periodic clus- 
tering of these events, as shown in chapter 5.2. 

Summarizing, the many abrupt climate changes 
characterising Earth's Ice Age epoch may reflect the 
interaction between Z and Earth during the last 3.2 
Myr. These interactions were global catastrophes, 
the worst of which resulted in a polar wandering. 
These particular events even irreversibly changed 
Earth's climate. The last polar shift event occurred 
only 11.5 kyr ago. As pointed out in chapter 2, 
the violence of this event is well documented by 
the observed sudden death of many animals, the ex- 
tinction of some of the larger mammal species and 
the sudden rise of the sea level Q. Possibly the 
traditions in old cultures of catastrophic events re- 
fer to the events mentioned above. Presently, we 
are enjoying a perfectly Milankovitch-driven mod- 
erate climate without the extreme and fast temper- 
ature variations observed prior to the last polar shift 
event. This may be due to the fact that planet Z 
is no longer influencing the Earth. We are confi- 
dent that this new and pleasant climate state, called 
Holocene, shall last for an indefinite time, provided, 
we are willing to stop our own "Greenhouse attack" 
on Earth's climate. 

Although we have presented a consistent expla- 
nation of the complex structure of Earth's climate 
variation during the last few Million years, we are 
aware that this rather exotic story has opened up 
more questions than it answered, because many of 
our arguments are based on considerations which 
are in terms of order of magnitude and very ele- 
mentary. Nevertheless, we believe that this proposal 
merits further investigation, in view of its basic sim- 
plicity. In particular, we encourage state of the art 
calculations which could narrow the allowed param- 
eter range of the hypothetical planet, or even dis- 
prove its existence. Of importance would be model 
calculations for the global climate for both pole po- 
sition, with and without the suggested additional 
Greenhouse effect. Our static approach to Earth's 
deformation during a fly-by should be replaced by 
a dynamical calculation. The actual amount of ma- 
terial that evaporates from a large body during its 
closest approach to the sun deserves a more detailed 
investigation. The sequence and efficiency of reac- 
tions initiated by an energetic atom, mostly Oxygen 
and Silicon, impinging on and stopped in the up- 
permost atmosphere, has to be further clarified, if 
necessary, with the help of appropriate experiments. 
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A determination of the maximum age of the ice at 
the Siberian side of the Arctic and/or at the At- 
lantic side of the Antarctic could provide a test of 
the polar shift theory. Planet Z must have left its 
chemical signature in polar ice. The material was 
deposited in atomic or molecular form, in contrast 
to dust which is mainly of terresrial origin. Sam- 
ples from the Moon could also reveal the presence 
of Z. These chemical traces might decide in which 
way Z disappeared from the solar system. If the 
dust cloud which envelops the inner solar system 
near the ecliptic plane originated from the last frag- 
ments of Z, then the intensity of the zodiacal light 
might show a measurable decrease in time pc| . We 
are optimistic, that the combined results of such 
studies will finally converge to a new understand- 
ing of Earth's Ice Age epoch. A negative answer 
implies that the strange variability of the climate 
during this epoch, as well as the stepwise cooling of 
Earth's climate beginning about 3.2 Mio years and 
its abrupt ending only 11.5 kyr ago remain unex- 
plained. We are sceptical that these problems can 
be solved by means of any of the presently favoured 
Milankovich-based theories alone. 
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